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If a body in a stream is provided with small ridges aligned in the local flow direction, a remarkable 
drag reduction can be reached under turbulent flow conditions. This surprising phenomenon is 
called the 'shark skin effect'. We demonstrate, that a reduction of resistance can also be reached in 
creeping flows if the ridges are aligned perpendicular to the flow direction. We especially consider 
in gravity-driven film flows the effect of the bottom topography on the mean transport velocity. 

PACS numbers: 47.15.Gf 



INTRODUCTION 

It is yet controversely discussed that the dermal surface 
morphology of sharks is in order to improve the sharks' 
swimming performance 0. Nevertheless, it is widely ac- 
cepted that for bodies in turbulent flows a reduction of 
skin friction by some percent can be reached if the surface 
of the body is provided with small ridges aligned in the 
local flow direction [^|^. This rather counter-intuitional 
phenomenon occurs in turbulent flows if the riblct spac- 
ing is smaller than the typical diameter of the streamwise 
vortices such that the vortices are forced to stay above 
the riblcts and touch only the tips of the riblcts. 

The criterion for turbulent flow conditions is based on 
the Reynolds number 



Re := 



gUL 



(1) 



which gives a comparison between inertial forces and fric- 
tion forces. By U and L a characteristic flow velocity and 
a characteristic length of the system are denoted, g and 
77 are the mass density and the dynamic viscosity of the 
fluid. At high Reynolds numbers, typically 10"* or larger, 
the flow is turbulent and vortices are created due to in- 
ertia. At small Reynolds number the flow is laminar and 
at vanishingly small Reynolds numbers Re <C 1 the flow 
is creeping. Hence, in creeping flows a 'shark skin effect' 
is not expected since inertia-induced vortices are absent. 

If the the riblets are aligned perpendicular to the flow 
direction, however, kinematically induced vortices can be 
created even in creeping flows where inertia is absent 0, 
1^. We especially consider a steady, gravity-driven fllm 
flow of an incompressible Newtonian fluid on an inclined 
plane. The bottom of the inclined plane is provided with 
periodic corrugations according to FIC^ The mean film 
thickness is denoted by H, the periodic length by A and 
the mean inclination angle by a. A Cartesian coordinate 
system is used with the x-axis placed at the mean level 
of the bottom contour, the y-axis in line with the ridges 
and the z-axis normal to the mean level of the bottom. 




FIG. 1: Sketch of the flow geometry 

MATHEMATICAL FORMULATION 
Basic equations and solving procedure 

The basic field equations are the continuity equation 
and the Navier-Stokes equations 



V • V 
- + (..V). 



0, (2) 



-\'p + riV^v + gg . (3) 



By g the gravity acceleration is denoted. The bottom 
contour is in general characterised by a 27r-periodic func- 
tion b{x) as z = 6(27ra-/A). Along the bottom the flow 
fulfills the no-slip condition 







(4) 



Assuming that the film height exceeds the wavelength 
of the bottom, the curvature of the free surface can be 
neglected Then, the free surface is given a& z = H. 



2 



At the free surface the dynamic boundary condition 

ez{[p-Ps]i-?/[V0i;+(V®i/)^]} = 0, (5) 

which is the equiUbrium between pressure jump at the 
surface and viscous forces, has to be fulfilled. By ps 
the pressure of the surrounding is denoted. If the spa- 
tial extensions of the bottom in y-direction are suffi- 
ciently large, side wall effects can be neglected and a 
two-dimensional flow geometry can be assumed. This 
allows for a representation of the velocity field in terms 
of a stream function, i.e. 



dip ^ dip _ 



dz 



dx 



(6) 



by which the continuity equation (|2Jl is identically ful- 
filled. Since due to the creeping condition Re ^ 1 the 
inertia terms at left hand of the Navier-Stokes equations 
are negligible, EQ. Q simplifies to the Stokes equations 

= -Vp + 77V^u + £»5 . (7) 

the solution of Stokes equations 



Reconsidering EQ 
is explicitly given as 



iP 



ggH^ sin a 



P = Ps 



27, 

- Q9 



z 



■nH^ sin a , , , 

{H - z) cos a + 

A 



(8) 
(9) 



where q{S^) and r(^) are holomorphic functions of the 
complex variable 



z -\-\x 
X 



(10) 



The symbols 5R and $5 denote the real and imaginary 
part of a complex expression, the prime derivation with 
respect to ^. By inserting the general solution (jH) EJ) in 
the boundary conditions |0][SJ), a set of equations for the 
boundary values of the two holomorphic functions q{S,) 
and r(^) is derived. Due to the periodicity of the flow a 
representation of the boundary values of g(^) and r(^) by 
Fourier series can be applied, which leads to an algebraic 
set of equations for the series coefficients. After trun- 
cating to a finite number of Fourier modes, its solution 
is determined by using computer algebra, e.g. Maple. 
The above method and the solving procedure is described 
in detail in j^. For the calculations presented here the 
series have been truncated at the mode number where 
the values of the coefficients fall below 10"^". This leads 
to truncation orders between 12 and 36. 
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FIG. 2: The three diflferent bottom shapes b±, 63 and 61 



narrow peaks. We approach this strongly idealised shape 
by three different trigonometric polynomials as shown in 
FIG. 121 At first, we consider a harmonic contour, i.e. 



-Acosx 



(11) 



with amplitude A as a comparatively smooth shape. In 
contrast to this, by the third order trigonometric poly- 
nomial 



63 (i) = - — [15 cos £ 
lb 



6 cos 2x + cos 32:1 



(12) 



an array of peaks is given. Since the derivatives of 63 (i) 
up to the 5th order vanish, a flat region between the 
peaks becomes apparent. The parameter A is again the 
amplitude, in the sense that fo3(7r) — ^3(0) = 2yl. The 
third shape shown in FIG. El which is defined as 



fcio(i) = 



A 



167960 cosx - 125970 cos 2i 

262144 L 

-1-77520 cos 3i - 38760 cos 4i + 1 5504 cos 5i - 4845 cos 6x 
+1140cos7i - 190cos8i + 20cos9i - coslOi , (13) 

is an array of significantly narrower peaks than in the 
case b3{x). Note, that the derivatives of bio{x) vanish 
up to the 19th derivative. Therefore, bio{x) is a good 
approximation for the idealised array of brush-like peaks 
by a continuous bottom shape. 

RESULTS 
Streamlines and vortex creation 



Bottom shapes 

For the modelling of wall roughness Panton sug- 
gests a 'brush model', i.e. an infinite array of equidistant 



Streamline patterns have been calculated for various 
shapes, amplitudes and film heights. For small ampli- 
tudes the streamlines follow the bottom contour, whereas 
fiow separation is observed if the amplitude exceeds a 
critical limit. As representative examples in FIG. the 



where the two-dimensional flow rate is given as 



(b) 

FIG. 3: Near-bottom streamlines of film flows over tfie 
sfiape biQ with amplitude (a) A = 0.17A and (b) A = 0.5A. 



streamhnes in the vicinity of the bottom are presented 
for the shape bio with two different amplitudes, namely 
A = 0.17 X and A = 0.5A. The film heights, H = 1.70A in 
(a) and H = 2.24A in (b), are chosen in order to receive 
the same flow rate in both cases. 

The critical amplitude for the primary flow separa- 
tion is A w 0.107A. Thus, in FIG. ISt flow separation 
is already apparent: A vortex pair has been created at 
the positions of maximum curvature. With increasing 
amplitude the vortices are growing, which leads to the 
merging of the two vortices to a single one. Such a case 
with a large single vortex which covers the major part 
of the region between two neighbouring peaks is shown 
in FIG. IHId. In this example we especially see a slightly 
curved separatrix passing nearly from tip to tip. This 
feeds the hope of a probable resistance reduction, since 
the vortices are supposed to act like 'fluid roller bearings'. 
By increasing the amplitude further, a secondary vortex 
pair is created at the critical amplitude A « 0.606A for 
secondary flow separation. 



Mean transport velocity 

We dcflne the mean transport velocity as 



H 

V= j udz = ^P{x,H)-'ilj{x,b{2nx/X)) . (15) 

b{2TTx/\) 

The quantity Ht, subsequently called the mean trans- 
port thickness, has to be understood as mean thickness 
of the part of the flow which contributes to the material 
transport, i.e. the film above the separation areas. If 
the separatrix of the primary vortex is given as z = s{x), 
xi < X < X2, the mean transport thickness results in 

Ht:^H~jJ [s{x) - 6(27rx/A)] dx . (16) 

The quantity Ht has to be carefully distinguished from 
the mean geometrical film thickness H of the entire film 
which includes the vortices. Furthermore, by 

n..4^'f (17, 

the reference thickness of a film fiow on a plane bottom 
with the same flow rate as the flow over the topography 
is defined 0. The three quantities TJq, H and Ht are 
illustrated in FIG. 01 Thus, for a fixed flow rate V , the 
comparison of the mean transport thickness Ht with the 
reference thickness Hq delivers an adequate measure for 
enhancement or reduction of the mean transport velocity 
in the fllm: In case of Ht > Hq, the mean transport veloc- 
ity is reduced, whereas Ht < Hq indicates enhancement 

of Ut. 

In FIG. |S1 the relative film elevation {Ht — Ho)/Ho 
is plotted versus the amplitude A for the three differ- 
ent bottom contours hi, 63 and 610. This parameter 
study has been carried out with a fixed flow rate of 
V = 9ggX^ sin a/ (Sry), corresponding to a reference thick- 
ness of Hq = 3A/2. Additionally, the onset of primary 
and further flow separation is indicated in the diagram. 
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FIG. 5: Relative film elevation vs. amplitude of flows over 
different shapes. 

From the beginning up to an amplitude of about A/6, 
the film elevation is monotonously increasing for all three 
different shapes, which indicates a reduction of the mean 
transport velocity due to the bottom corrugations. Note, 
that within this parameter regime no positive effect can 
be expected since no vortices are present. However, the 
curves reach maxima slightly after the primary vortex 
generation and pass then into a monotonous decrease 
due to reduction of friction by vortices, which act like 
fluid roller bearings. Obviously, both the height of the 
maximum as well as the decrease after the vortex genera- 
tion are more pronounced for bottom shapes with sharper 
peaks. For the curve associated to the bottom 6io the 
film elevation becomes negative for high amplitudes, i.e. 
the mean transport velocity exceeds the corresponding 
mean transport velocity of the flow over a plane bottom. 
At the highest amplitude considered in our calculations, 
A = 0.9A, this enhancement of the mean transport ve- 
locity reaches 0.88%. 

The increase of the mean transport velocity indicates 
an improved mass transport in the volume. In contrast 
to this, the surface velocity remains nearly unchanged. 

CONCLUSION 

The parameter studies on the three different shapes 
&i, &3 and &10 revealed a noticeable effect of the bot- 



tom topography on the material transport in a creep- 
ing film flow: For the shapes hi and 63 we found a de- 
crease of the mean transport velocity compared to the 
flow over a plane bottom, whereas for the shape &10 with 
the sharpest peaks an increase of the mean transport ve- 
locity becomes apparent at sufficiently high amplitude. 
Thus, a comparison to the 'shark skin effect', which has 
been successfully applied to ships and airplanes for drag 
reduction, is near at hand. The present effect, however, 
is essentially different from the popular shark skin effect: 
The increase of the mean transport velocity has been cal- 
culated for creeping flows, whereas the shark skin effect 
occurs in turbulent flow at Reynolds numbers w 10^-10^. 
Furthermore, the rippled structures of the shark skin are 
directed longitudinal to the flow, not transversal as it is 
the case here. Finally, the responsible mechanism for the 
shark skin effect is according to Q and Q the control 
of the streamwise vortices in the turbulent flow, which 
has been created by inertia. In contrast to this, the rip- 
pled bottom structure in creeping film flows enforce the 
creation of vortices which act on the flow like a kind of 
'fluid roller bearing'. Nevertheless, a common feature of 
the shark skin effect and the effect observed in the present 
paper is the reduction of resistance in the flow by means 
of rippled wall structures. 

The highest relative increase of the mean transport 
velocity is in our calculations 0.88%. This is less than 
the shark skin effect in turbulent flow, which is supposed 
to be at most 10% 0- It is an open question up to which 
extend the value of 0.88% can be improved by varying 
the relevant parameters, especially the bottom shape. 
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